A New Search Paradigm for Correlated Neutrino 
Emission from Discrete GRBs using Antarctic 
Cherenkov Telescopes in the Swift Era 

v©"- Michael Stamatikos for the IceCube 2 Collaboration* and David L. Band f 

O 

' * Department of Physics, University of Wisconsin, Madison, WI 53706, USA^ 

^NASA/GSFC, Greenbelt, MD 20771, USA and JCA/UMBC, Baltimore, MB 21250, USA 4 

<D 

Ph ' Abstract. We describe the theoretical modeling and analysis techniques associated with a prelim- 

, inary search for correlated neutrino emission from GRB980703a, which triggered the Burst and 

^vq ■ Transient Source Experiment (BATSE GRB trigger 6891), using archived data from the Antarctic 

Muon and Neutrino Detector Array (AMANDA-BIO). Under the assumption of associated hadronic 
acceleration, the expected observed neutrino energy flux is directly derived, based upon con- 
fronting the fireball phenomenology with the discrete set of observed electromagnetic parameters 
of GRB980703a, gleaned from ground-based and satellite observations, for four models, corrected 
for oscillations. Models 1 and 2, based upon spectral analysis featuring a prompt photon energy fit 
to the Band function, utilize an observed spectroscopic redshift, for isotropic and anisotropic emis- 
sion geometry, respectively. Model 3 is based upon averaged burst parameters, assuming isotropic 
emission. Model 4, based upon a Band fit, features an estimated redshift from the lag-luminosity 
relation with isotropic emission. Consistent with our AMANDA-II analysis of GRB030329, which 
resulted in a flux upper limit of ~ 0.150GeV/cm 2 /s for model 1, we find differences in excess of an 
■ order of magnitude in the response of AMANDA-B 10, among the various models for GRB980703a. 

Implications for future searches in the era of Swift and IceCube are discussed. 
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§ : 1. GRB030329 & THE CASE FOR A NEW PARADIGM 
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Canonical fireball phenomenology, in the context of hadronic acceleration, predicts cor- 
related MeV to EeV neutrinos from gamma-ray bursts (GRBs). Ideal for detection are 
~ TeV-PeV muon neutrinos, which arise as the leptonic decay products of photome- 
son interactions (p + y —> A + — > 7t + + [n] — > /i + + v M — > e + + v e + v M + v M ) within 
the internal shocks of the relativistic fireball. Since the prompt y-rays act as the am- 
bient photon target field, these neutrinos are expected to be in spatial and tempo- 
ral coincidence, with inverted 5 energy spectra (see Equation El, which trace the pho- 
ton energy spectra (see Equation [B, due to the intrinsic threshold requirement that 
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5 Consequently, e v °= £y l « £,,, as regulated by the normalization (A V/1 ) and proton energy efficiency {f K )- 



4(1 +z) 2 £ p £y > — inV) r^ ulk . Constraints imposed by coincidence are tantamount 
to nearly background-free searches in neutrino observatories such as AMANDA and 
IceCube (see Table[T])- The former, which has been calibrated with atmospheric neutri- 
nos, has demonstrated the viability of high energy neutrino astronomy using the ice at 
the geographic South Pole as a Cherenkov medium since 1997. The latter, AMANDA' s 
km-scale successor, is currently under construction with anticipated completion by 2010. 

A positive detection of such high energy neutrinos would confirm hadronic accel- 
eration in the relativistic GRB-wind, providing critical insight to the associated micro- 
physics of the fireball, while possibly revealing an astrophysical acceleration mechanism 
for the highest energy cosmic rays. AMANDA analyses have used a diffuse formulation 
HI], predicated on an ensemble of average GRBs, to produce the most stringent upper 
limits upon correlated multi-flavored neutrino emission However, over 30 years of 
ground-based and satellite electromagnetic observations have documented the following 
facts: (i) the electromagnetic parameters of GRBs are characterized by distributions, of- 
ten spanning multiple orders of magnitude, which differ both among and within bursts of 
different classes (i.e. short, long, x-ray rich, etc.), whose deviation from averaged values 
is often not accommodated by the inherent uncertainty of measurement; (ii) GRB satel- 
lite detectors exhibit a large dynamic range of thresholds and sensitivities which impart 
statistical sample bias via selection effects. These facts render the notion of an average 
GRB incompatible with the observational record. Furthermore, it has been argued that 
electromagnetic variations maylead to variations in the number of expected neutrino 
events associated with GRBs [3]. This reasoning has led to a new modeling paradigm 
for correlated neutrino emission searches, which is based upon the notion of a discrete 
GRB, i.e. one that is described by a unique set of electromagnetic parameters. 

The quantitative effects on the expected neutrino number, energy and subsequent con- 
straints upon astrophysical models (due to null detection), based upon multiple models 
of emission geometry and electromagnetic characterization, were initially illustrated in 
an analysis of GRB030329 with AMANDA-II 0]. For GRB030329, a peak effective 
area for muon neutrinos of ~ 80 m 2 and ~ 700 m 2 at ~ 2 PeV and an effective area for 
muons of ~ 100,000 m 2 and ~ 1 km 2 at ~ 200 TeV were achieved for AMANDA-II 
and predicted for IceCube, respectively. Principal results, including neutrino flux up- 
per limits for each model tested, are summarized in Table d Further details regarding 
this analysis may be found elsewhere 0]. Here, supplementary to our conference pre- 
sentation, we extend the paradigm by taking a first look at GRB980703a, one of the 
BATSE bursts currently under investigation using AMANDA archived data from 1997- 
2000 using models 1 (discrete-isotropic), 2 (discrete-jet), 3 (average-isotropic) 

and 4 (lag-isotropic). Model 4 demonstrates that the lag-luminosity relation [7] provides 
a reasonable redshift estimate. 



2. GRB980703A: EXTENDING THE DISCRETE GRB PARADIGM 

Since the expected neutrino emission relies on the microphysics associated with the 
prompt y-ray emission phase, a complete electromagnetic characterization of the GRB 
is compulsory. Spectral analysis, resulting in a prompt y-ray photon energy spectrum, 



TABLE 1. Summary of results for GRB030329 fl. 



Model n b , n' b * N s , n s , rif n ohs , n' obs ** Flux Upper Liming 

1 17.44,0.06 0.1308,0.0202,0.0156 15,0 0.150 

2 17.44,0.06 0.0691,0.0116,0.0092 15,0 0.039 

3 17.44,0.06 0.0038,0.0008,0.0006 15,0 0.035 



* Number of background events expected during a 40 second on-time search window before (nb) and 
after (n' b ) quality selection (optimized for discovery 0|), including restriction to a search bin radius 
(space angle between the reconstructed muon trajectory and the GRB's position) of 1 1.3°. 
^ Number of neutrino signal events expected on-time for IceCube (N s ) and AMANDA-II (n s , n' s ). 
** The number of observed events in AMANDA-II before (n^) and after (n' , ) quality selection. 

$ Based upon null detection in AMANDA-II. The effects of neutrino flavor oscillations have been 
included. For more details on GRB030329's electromagnetic and neutrino parameterization, see 0|. 
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FIGURE 1. Left Plate - Prompt photon energy spectral fit to the Band function (see Equation^ f° r 
GRB980703a using the observed discrete electromagnetic parameters (solid black - used in models 1, 2, 
and 4), and average values (solid red - used in model 3). Inset illustrates the regime of the photon break 
energy (e* ) for the discrete parameter fit. Note that the disagreement between the curves is not resolved by 
the uncertainty in the discrete fit, illustrated as 1 a dashed confidence bands (assuming propagated errors 
with trivial covariance). See Tables [2] and [3] for fit parameters. Right Plate - The prompt photon energy 
spectral Band fit, convolved with BATSE detector response (solid line), is compared to the measured 
BATSE photon count rate (data points with error bars) for the entire (fluence) emission of GRB980703a. 



illustrated in Figured was performed via convolving an assumed spectral model with 
the BATSE detector response and comparing the fit with observed data. Our spectral 
assumption is characterized by an empirical model known as the Band function (a]: 
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Multi-wavelength afterglow observations have resulted in a more comprehensive elec- 
tromagnetic characterization of GRB980703a, such as an observed spectroscopic red- 



shift, Zobsi via analysis of the host galaxy II 1011 . Since the blast wave evolution is sensitive 
to intrinsic burst properties such as the explosion energy, emission geometry and cir- 
cumstellar density, n, the afterglow has been used to deduce anisotropic emission from 
a jet break time, tj et , and the relative energy fractions imparted to both electrons, £ e , 
and the magnetic field, £g. The best fit solution is consistent with the canonical fireball 
model, with a jet emission geometry that is inferred by assuming a relativistic outflow, 
beamed into an ambient medium of constant density, within a collimated jet of half an- 
gle, 0jet JO- A beaming fraction correction, fs, is then used to correct isotropic values 

for luminosity, Ly°, and energy, Ey°, into beam-corrected (jet) values of Ly l and Ey , 
respectively. Electromagnetic observables are summarized in Table while the electro- 
magnetic parameterization of each model is given in Table 
TABLE 2. Observed electromagnetic properties of GRB980703a (BATSE Trigger #: 6891). 



Parameter(s) 


Value 


Reference 


RA, DEC, C7 R (°) [J2000] 


359.777775417,8.585303861, 1.38 x 10~ 7 


[12J 


T (UTC S ), T 90 (s)* [50-350 keV] 


15765.22, 411.65 ±9.27 


[13J 


F y Total (ergs/cm 2 ) [50-350 keV] 


(6.22 ±0.43) x 10~ 5 


[13] 


4>£ eak (ergs/cm 2 /s) [20-1000 keV] 1 " 


(1.56 ±0.07) x 10~ 6 


® 


a, jS, £° (keV) [50-350 keV]** 


-1.02 ± 0.20, -2.33 ± 0.57, 223 ± 96 


16] 


e\, e ! y (keV) [50-350 keV] 


293 ±185, 219 ±104 






0.9660 ±0.0006 


[10J 




6219±2843Mpc- (1.92 ±0.88) x 10 28 cm 


® 


£e, £b 


0.27 ± 0.03, 0.0018+JJ'^ 


U4] 


tjet (days) 


3.43±0.50 


[11] 


n (cm -3 ) 


28±10 


111] 


fife* 


- (10. 10 ±1. 36)° w (0.176 ± 0.024) rad 




f B = 1 - COS ©jet 


0.016 ±0.004 


Llil 



* The start of the Tgo interval, with respect to the trigger time (T), was -6.14 seconds 1 13]. 

f 0jet « 0.101 rad (^) 1 (&) 4 (_^) 4 [Op] "* (^) with § ~ 0.2018S (see Q). 
' We take the energy band pass of 20-1000 keV as bolometric. 

** A y = (1.97±0.43) x 1(T 3 photons/cm 2 /keV/s, xl ~ 0.80, and signal/noise ~ 4.11 (see Equation^. 

* Based upon an average from emission and absorption (Doppler redshift) lines from the host galaxy full . 

§ Acdm cosmology: H D = 72 ± 5 km/Mpc/s, £2 m = 0.29 ± 0.07, £1 A = 0.73 ± 0.07 O, is utilized throughout. 

The neutrino spectral parameterization is given in Equation 1 4]. Values from Ta- 
bles El and El have been substituted into Equation|2]to produce the values of the neutrino 
parameters and the number of neutrino events expected in AMANDA-B10 (via simula- 
tion) for models 1-4, given in Table0 The response of AMANDA-B10 to GRB980703a 
models 1-4, including effective areas for neutrinos and muons, is illustrated in FigureEl 
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TABLE 3. Electromagnetic parameterization for GRB980703a neutrino models. 



Parameter 



Model 1* Model 2 1 " Model 3* 



Model 4* 



z 




See Table]2\ 


See Table\2\ 


~ 1 


0.6378+S 


d L (Mpc) 




See Table\2\ 


See Table\2\ 


6491 ±2887 


3730+H 


L r (10 51 erj 


;s/s) [20-1000 keV] 


7.21 ±6.60 


0.11 ±0.11 


10.08 ±8.97 




E r (10 52 erj 


is) [20-2000 keV] 


13.77 ±12.68 


0.21 ±0.21 


1.51 ±1.35 


6 57+ 5 26 






See Tabled 


See Tabled 


33+ 67 

u -"-0.33 


33+ 67 


£b 




See Table\2\ 


See Table\2\ 


u - JJ -0.33 


11+ - 67 

U ""-0.33 



* Based upon discrete electromagnetic parameters, assuming isotropic emission. 
' Based upon discrete electromagnetic parameters, assuming beamed (jet) emission. 



Using averaged GRB values: a~ — 1, J3 ~ —2, £y ~ e£ ~ e£ 



1 MeV, FT° tal . 



- 6 x 10 ergs/cm 



2 „. 



and 4>^ eak ~ 2 x 10 6 ergs/cm 2 /s Q[. Under the assumption of isotropic emission. 
* Isotropic emission and estimated z via the lag-luminosity method, with T a = seconds Q]. 

TABLE 4. Prompt muon neutrino flux parameterization and event number for GRB980703a. 



Parameter* 



Model 1 Model 2 



Model 



Model 4 



Tbuik ^ 276 



^.52^2<mW(1+^) 



293 ±45 



146 ±23 



310 ±46 [300] 



243 



+33 
-17 



h - 0-2 x -r— 

i 2.5 f v- 
Fyfx 



v.-2 £ ^M,.v( 1 + ;: ) 



0.34 ±0.43 0.09 ±0.11 0.11 ±0.12 [0.20] 0.34+^ 



0.40 

28 



^sipife(10- 6 GeV/cm 2 /s) 



7 x io 5 r| 5 

fl +z) 2 



(10 s GeV) 



6.49±8.27 1.62±2.13 
5.29±3.71 1.32 ±0.93 



0.16 



-0.37 



£ t w i^gA^L (ifjs GeV) 



(l+z)v/£SA/iy-52 



5.37 



+4.15 



^ AMANDA-B10 (JO -4 )* 



-4.13 

6.08 



2.68 



-2.16 
-2.15 

2.79 



J -0.24 

1.68±0.50[1] 



0.4612-25 [0.1] 



5.27 
5.25 



J -0.47 

0.32 



0.40 



12.36 
6.81 

+3.88 
-3.51 

-0.64 



-0.32 
4.95 



Where L 7 = L y52 ■ 10 52 ergs/s, T = T 2 . 5 • 10 2 ' 5 , t v 



1 MeV, and e" v lax = e^ eV • 100 MeV. 



: t v — 2 • 10ms, £y — ^y.MeV ' 1 1VA ^ v ' c 7 

' Bracketed values indicate nominal values used in average GRB parameterization Q[. 

** Super-Eddington luminosity within a compact source requires a lower bound to ensure transparent optical depth. 

* Events expected in AMANDA-B10, on a background of ~8. 17, for an on-time search window of 420.92 seconds 
within a search bin radius of ~ 9.5°. Estimated event selection attenuates signal efficiency by ~ 75%, while rejecting 
~ 99.7% of the background. These values are preliminary estimates, a full analysis is currently in progress. 



3. DISCUSSION & FUTURE OUTLOOK IN THE SWIFT ERA 

Consistent with our results for GRB030329, we find that the most critical parameters, 
which translate into an observable variation in detector response, are the electromagnetic 
fluence, F r , and spectral characterization in the vicinity of the photon break energy, et. 
The former is related to the number of neutrinos expected in the detector, while the 
latter affects the mean neutrino energy of the events. These effects on neutrino energy 
flux, which exceed one and two orders of magnitude in the comparison of models 1 and 
3, in GRB980703a and GRB030329, respectively, may be directly traced back to the 
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FIGURE 2. Left Plate - Prompt neutrino energy flux {Upper panel) and AMANDA-BIO detector 
response (Lower panel), for GRB980703a models 1, 2, 3, and 4, indicated by solid black, dashed blue, 
dot-dashed red, and hatched magenta curves, respectively. Right Plate - Effective area of AMANDA-BIO 
for muon neutrinos (Upper panel) and muons (Lower panel - as a function of energy at closest approach 
to the detector). Color code is identical to left plate. See Table@]for selection and signal efficiency details. 



variance in their electromagnetic characterization, which directly affects the constraints 
placed upon astrophysical models in the case of null detection. The consequences of 
this reality are both unequivocal and apropos, since a hallmark of the Swift era is the 
acquisition of a more complete electromagnetic characterization of fewer bursts (relative 
to the age of BATSE). With an order of magnitude increase in effective area for GRBs 
0], realized by km-scale detectors such as IceCube, it is most likely that either the 
first evidence for correlated neutrino emission or the first real constraints on associated 
hadronic acceleration will come from the analysis of an exceptional (local) discrete 
GRB, rather than an aggregate of hundreds with average emission. It is interesting to 
note that model 4, based upon an estimated redshift from the lag-luminosity relation, 
was consistent with model 1, which was based upon the observed redshift. Future work 
includes the analysis of a subset of BATSE GRBs from 1997-2000 JUl, using Band 
function fits and (lag-luminosity relation) estimated redshifts. Ultimately, a synergy of 
gamma-ray and neutrino astronomy may be realized, within the context of GRBs, via 
multi-wavelength and multi-messenger correlative observational campaigns in an era of 
superior scientific instruments such as IceCube and Swift. 



REFERENCES 

1 . E. Waxman, and J. Bahcall, Phys. Rev. D 59, 023002 (1999). 

2. A. Achterberg, et al., astro-ph/0509330 (See also K. Kuehn et al., these proceedings.). 

3. D. Guetta, et al., Astroparticle Physics 20, 429-455 (2004). 

4. M. Stamatikos, J. Kurtzweil, and M. Clarke, preprint (astro-ph/05 10336) (2005). 



5. M. Stamatikos, et al., "Gamma-Ray Bursts: 30 Years of Discovery," AIPC 727, 2004, pp. 146-149. 

6. M. Stamatikos, D. Band, D. Hooper, and F. Halzen, In preparation, to be submitted to ApJ. (2006). 

7. D. Band, et al., ApJ 613, 484^91 (2004). 

8. G. Hill, J. Hodges, B. Hughey, A. Karle, and M. Stamatikos, "PHYSTAT," Oxford, 2005. 

9. D. Band, et al., ApJ 413, 281-292 (1993). 

10. S. Djorgovski, et al., ApJ 508, L17-L20 (1998). 

11. D. Frail, et al., ApJ 590, 992 (2003). 

12. G. Taylor, et al., GCN GRB Observation Report 152 (1998). 

13. (2003), URL |http : / / www .batse.msfc.nasa. gov/ bat se/grb/ catalog/ cur rent /| 

14. S. Yost, et al., ApJ 597, 459-473 (2003). 

15. A. Friedman, and J. Bloom, ApJ 627, 1-25 (2005). 

16. D. Spergel, et al., ApJ Supp. Ser. 148, 174-194 (2003). 



